MATH 18.152 COURSE NOTES - CLASS MEETING # 2 



18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 2: The Diffusion (aka Heat) Equation 



1. Introduction to the Heat Equation 
The heat equation for a function u(t,x), x d = (x 1 , ••■,x n ) e IR n , is 

(1.0.1) u t -DAu = f(t,x). 

Here, the constant D > is the diffusion coefficient, f(t,x) is an inhomogeneous term, and A is the 
Laplacian operator, which takes the following form in Cartesian coordinates: 

n 

(1.0.2) A d = cf £df. 

i=i 



Equation (1.0.1) is first-order and linear. 



2. A SIMPLE MODEL OF HEAT FLOW THAT LEADS TO THE HEAT EQUATION 

We now give an example of a simple model of heat flow that leads to the heat equation. Consider 
a homogeneous, isotropic solid body Scl" (n = 3 is the physically relevant case) described by the 
following physical properties: 

(2.0.3) p d = mass density ~ [mass] x [Volume] -1 = constant, 

(2.0.4) e(t,x) = thermal energy per unit mass ~ [energy] x [mass] . 

Let's also assume that heat is supplied to the body by an external source which pumps in heat at 
the following rate per unit mass: 

(2.0.5) 8% ~ [energy] x [time]" 1 x [mass] -1 . 

The total thermal E(t; V) energy contained in a body sub-volume V c B at time t is the integral 
of e(t, x) over V ■ 

(2.0.6) E(t;V) d = f pe(t,x)d n x. 

Jv 

The rate of change of the total energy contained in V is 

(2.0.7) ^E(t; v ) = ^ t f v P e (t, x) d n x = pd t e(t, x) d n x. 



In (2.0.7), we have assumed that you can differentiate under the integral; we can do this when 
e(t, x) is a "nice" function. We will be more precise about the meaning of "nice" later in the course. 

l 
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Let's now address the factors that can cause 4iE(t; V) to be non-zero. That is, let's account for 
the factors that cause the energy within the volume V to change. In our simple model, we will 
account for only two factors. First, by integrating ( |2.0.5 ) over V, we deduce the rate of energy 
pumped into the sub-volume V by the external source: 

(2.0.8) / p&(t,x) d n x ~ [energy] x [time]" 1 . 

Jv 

Second, we will also assume that heat energy is flowing throughout the body, and that flow can 
be modeled by a heat flux vector q 

(2.0.9) q ~ [energy] x [time]" 1 x [area]" 1 , 

which specifies the direction and magnitude of heat flow across a unit area. That is, if da c dV is 
a small surface area with outward unit-normal N, then q • N is the energy flowing out of the small 
surface. Thus, the rate of heat energy flowing into V is 

(2.0.10) -/ q-Nda = - V ■ qd n x ~ [energy] x [time]" 1 , 

Jdv Jv 

where the equality follows from the divergence theorem. 

We will connect the various energies together by assuming the following energy conservation 
"law:" The rate of change of total energy in the sub-volume V is equal to the rate of heat energy 



flowing into V + rate of heat energy supplied by the external source. Using (2.0.7), (2.0.8), and 



(2.0.10), we see that this "law" takes the following form in terms of integrals: 



(2.0.11) ^ pd t e(t, x) d n x = - ^ V • q d n x + ^ p$£ d n x. 

Since the above relations are assumed to hold for all body sub-volumes V, the integrands must be 
equal (again, as long as they are nice): 

(2.0.12) pd t e(t,x) = -V-q + p^. 



2.1. Fourier's law. In order to turn (2.0.12) into a PDE that we can study, we need to make an- 
other assumption about e(t, x), q, and their relation to the temperature u(t, x). Fourier hypothesized 
the following "Fourier's Law of heat conduction:" 

(2.1.1) q(t,x) = -KVu(t,x), 

where k > is the thermal conductivity, and Vw d = f (diu,---,d n u) is the spatial derivative gradient 
of the temperature u(t,x). We will assume that k is a constant. Recall that at each fixed t, 
Vu(t,x) points in the direction of maximal increase and that S/u(t,x) is perpendicular to the level 



sets {x | u(x) = constant}. Thus, (2.1.1) states that heat flows "from hot to cold" (i.e. towards 



decreasing temperature) and that the flow is perpendicular to the surfaces of constant temperature. 



Remark 2.1.1. (2.1.1) is NOT A FUNDAMENTAL LAW OF NATURE] It is a simple but rea- 



sonable (under certain circumstances) model! 
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We need one more assumption in order to derive our PDE - we need to relate e(t, x) to u(t, x). We 
will assume a very simple model, which is experimentally verified by many substances in moderate 
temperature ranges: 

(2.1.2) e = c v u. 

Here, c v > is the specific heat at constant volume. We also assume that c v is constant. Like 



many of our previous assumptions, (2.1.2) is also just a simple model, and not a fundamental law 
of nature. 



Finally, we combine (2.0.12), (2.1.1), and (2.1.2), and use the identity V ■ = Am, thus arriving 



at 



(2.1.3) 



d t u(t,x) 



K 1 

— Au + — < 



c v p 



This is the heat equation (1.0.1) with D = ^- and / 



i 



3. Well-posedness 

Remember, one of the main goals of PDE theory is to figure out which kind of data lead to a 
unique solution. It is not always obvious which kind of data we are allowed to specify in order to 
solve the equation. When we have a PDE and a notion of data such that the data always lead to a 
unique solution, and the solution depends "continuously" on the data, we say that the problem is 
well-posed. 

3.1. Dirichlet boundary conditions. Let's study Dirichlet boundary conditions for the heat 
equation in n = 1 dimensions. Think of a one-dimensional rod with endpoints at x = and x = L. 
Let's set most of the constants equal to 1 for simplicity, and assume that there is no external source 
pumping energy into the rod, i.e., that there is no inhomogeneous term /. 

Then we could, for example, prescribe the temperature of the rod at t = (sometimes called 
Cauchy data) and also at the boundaries x = and x = L for all times t € [0, T] : 



' d t u - Dd 2 x u = 0, (t, x) e (0, T) x (0, L), 
(3.1.1) I u(0,x)=g(x), xe[0,L], (Cauchy data), 

u(t,0) = ho(t), u(t,L) = h L (t), t>0, (Dirichlet data). 

As we will see, under suitable assumptions on the functions, g,h ,h L , these conditions lead to a 
well-posed problem. 

3.2. Neumann (N for Normal!) boundary conditions. Instead of prescribing the temperature 
at the boundaries, let's instead prescribe the inward rate of heat flow (given by Fourier's law with 
k = 1) at the boundaries: 



' d t u-Dd%u = 0, (t,x) e (0,T) x (0,L), 
(3.2.1) I u(0, x) - g(x), (Cauchy data), 

-d x u(t,0) = h (t), d x u(t,L) = h L (t), (Neumann data). 

Under suitable assumptions on the functions, g,h ,hL, these conditions also lead to a well-posed 
problem. 



4 
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3.3. Robin boundary conditions. We can also take some linear combinations of the Dirichlet 
and Neumann conditions: 



(3.3.1) 



' d t u-Dd%u = 0, (t,x) e (0,T) x (0,L), 
u(0,x) = g(x), (Cauchy data), 

-d x u(t, 0) + ocu(t, 0) = ho(t), d x u(t, L) + <xu(t, L) = h L (t), 



(Robin data), 



where a > is a positive constant. Under suitable assumptions on the functions, g,ho,hi, these 
conditions also lead to a well-posed problem. 



3.4. Mixed boundary conditions. The above three boundary conditions are called homogeneous 
because they are of the same type at each end. It is also possible to prescribe one condition at 
one endpoint, and a different condition at the other endpoint. These are called mixed boundary 
conditions. These conditions also lead to a well-posed problem. 



4. Separation of variables 

We now discuss a technique, known as separation of variables, that can be used to explicitly 
solve certain PDEs. It is especially useful in the study of linear PDEs. Although this technique is 
applicable to some important PDEs, it is unfortunately far from universally applicable. 

In a nutshell, the separation of variables technique can be summarized as: 

• Look for a solution of the form u(t,x) = v(t)w(x). 

• Plug this guess into the PDE and hope that the PDE forces the functions v and w to be 
solutions to ODEs that can be solved without too much trouble. 

As we will see, when one tries to apply this technique, one quickly runs into difficulties that are 
best addressed using techniques from Fourier analysis. We don't have time right now to give a 
detailed introduction to Fourier analysis, but we will return to it later in the course if time permits; 
at the moment, we will only show how to use some of these techniques, without fully justifying 
them. 

A great way to illustrate separation of variables is through an example. Let's try to solve the 
heat equation problem with homogeneous (i.e., vanishing) Dirichlet conditions 



(4.0.1) 



' iH - u xx = 0, (t,x) € (0,T] x [0,1], 
w(0, x) = x, x e [0, 1], 
u(t,0) = 0, u(t,l) = 0, 



by separation of variables. 

Remark 4.0.1. Note that such a solution cannot possibly be continuous at the point (0, 1). 



We plug in the form u(t,x) = v(t)w(x) into (4.0.1) and discover that 



(4.0.2) 



v'(t) _ w"(x) 
v(t) w(x) 
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This should hold for all t, x. It therefore must be the case that both sides are equal to a constant, 
which we will call A. We then have 



(4.0.3a) v'(t)=Xv(t), 
(4.0.3b) w"{x)=\w{x). 

Furthermore, w(0) = w(l) = by the boundary conditions. 

Let's address v first, since it requires less work to deal with than w. If A € K, then (4.0.3a) can 
be generally solved: 

(4.0.4) v(t) = Ae M 

for some Ael. 

In contrast, the study of w(x) splits into three cases: 

• A = 0. Then w(x) = Bx + C for some 5,C e 1. The boundary conditions imply that C = 
and B + C = 0, so that B = C = 0. Thus, this solution is not very interesting. 

• A > 0. Then w(x) = Be s ^ Kx + Ce~^ x for some 5,CeM. The boundary conditions imply that 
B + C = 0, and Be^ + Ce~~^ = 0, which forces B = C = 0. This solution is also not very 
interesting. 

• A < 0. Then w(x) = B sm(^J\A\x) + C cos(a v /|a|x) for some B,Cel. The boundary condition 
w(0) = forces C = 0, so w(x) = Bsin(Vhx). The boundary condition w(l) = then forces 

A = -7r 2 m 2 for some m € Z + , where Z + d = f the set of non-negative integers. The A are 
called eigenvalues, and the corresponding w m are the corresponding eigenvectors. Equation 



(|4.0.3al) is called an eigenvalue problem corresponding to the linear operator C d = <9 2 



We have shown that the only solutions w are of the form w m (x) = Bsm(27rmx), m € Z + . Using also 



(4.0.4) and the fact that A = -7r 2 m 2 for our solutions, we have produced a family of solutions to the 



heat equation d t u - <9 2 w = that satisfying the boundary conditions: 
(4.0.5) u m (t, x) = e - mVi sin(mTrx), A m e K, m e Z + . 

But we haven't yet satisfied the initial condition u(0,x) = x. To do this, we could try using the 
superposition principle: 

oo 

(4.0.6) u(t,x) = Y, AnU m (t,x). 

m=l 

We would have to solve for the A m to achieve the desired initial condition u(0,x) = x. 

Here is a list of things we would have to do to fully solve this problem using this technique: 

(1) Find plausible A m . 



(2) Show that the infinite sum (4.0.6) converges. 

(3) Show that the infinite sum solves the heat equation. 

(4) Show that u(t,x) satisfies the boundary conditions. 

(5) Check that lim t ^o+ u(t, x) = u(0, x) = x. We also have to investigate in which sense this limit 
may or may not hold. We already know that this equality cannot hold pointwise at the 
point (0, 1). 
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(6) Show that there can be no other solution with these initial/boundary conditions (unique- 
ness). 



Let's deal with (1) first. If (4.0.6) holds, then at t = : 



(4.0.7) x = u(0, x) = Y, A m u m (0, x)= Y, A m sm(mnx). 



m=l m=l 



This is a Fourier series expansion for the function f(x) = x on the interval [0, 1]. 

It is helpful to think of a function f(x) as a vector in an infinite dimensional vector space and the 
sin(m7rx) as basis vectors (however, it is not trivial to show that they form a basis...). Furthermore, 
if we introduce the dot product 

(4.0.8) {f(x),g(x)) d i f f f(x)g(x)dx, 

J[0,1] 

then the basis vectors are orthogonal (do the computation yourself!): 

i a r, r,\ / • / \ • / \ \ I 1/2 if m = n 

(4.0.9) (sm{mnx),sm{7rnx)) = { „'. r 

This suggests that the following heuristic computations might be able to be made completely 
rigorous: 

p oo 

(4.0.10) / f(x) sin(7mx) dx = (f(x),sm(imx)) = (V A m sin(m7nr), sin(7rnx)) 

^[0,1] ' m= i 

oo 

= Y] (A m sin(m7rx), sin(7rnx)) 

m=l 

= -A n . 

Applying this to our function f(x) = x, we integrate by parts to compute that 

r 2 _ 2 r 2 

(4.0.11) A m = 2 / xsm(miTx) dx = x cos(mirx)\^lQ + / cos(irnx) dx = (-l) m+1 . 

J [o,i] mn mix J [o,i~\ mn 

We now hope that our solution is: 

oo o 

(4.0.12) u(t, x)=Y (-i)m+i e -m 2 nH S i n ( m7rx ). 

Remark 4.0.2. The individual terms (-l) m+1 e~ m27T2t — sm(mirx) are sometimes called the modes 

v ' mix v / 

of the solution. Note that each mode is rapidly decaying at an exponential rate as t -> oo. Further- 
more, the infinite sum Em=i( _i ) m+le "^"^r^ sin(m7ra:) also decays exponentially in time. Later 
in the course we will study the heat equation on all of R, and we will once again see that under 
suitable assumptions, solutions to the heat equation tend to exponentially decay in time. However, 



if we had non-zero Dirichlet conditions for the problem (4.0.1), then the solution might not decay 
to 0, but instead to some other state. 
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Let's now answer some of the remaining questions from above. 



(2) Thanks to the rapidly decaying in m factor e~ m 71 *, for any t > 0, the series (4.0.12) can be 
seen to uniformly converge for x e [0, 1] using one of the standard convergence arguments from 
analysis (carefully work through this argument yourself; pg. 9 of your book might be a helpful 



reference). The argument for t = is much more subtle and is addressed in Theorem 4.1 below. 



(3) We already know that each mode in (4.0.12) solves the heat equation. So what about the 



infinite sum? Again, for any t > 0, the e~ m27T 1 factor plus standard results from analysis allow us 
to repeatedly differentiate the series term- by-term in both t and x (work through this yourself). 
In particular, the series is smooth (i.e., infinitely differentiable in all variables) for any t > 0. In 
particular, for t > 0, we have that 

oo O oo 

(4.0.13a) d t u= Y, d t [(-l) m+1 e- m2n2t sin(mTrx)] = Y (-l) m rrme- m2 * 2t sin(m7rx), 

m=l m7T m=l 

oo o oo o 

(4.0.13b) dlu = T d 2 x [(-l) m+1 e- m2 * 2t sin(m7ra;)] = Y {-l) m m^e- m ^ 2t sin(m7ra:), 

rT^l ™ r^l m7T 

which shows that -d t u + d^u = 0. 

(4) The fact that u verifies the correct Dirichlet conditions at x = and x = 1 follows from the 
fact that each of the modes does. 

The remaining two questions require more work. We first quote the following theorem from 
Fourier analysis to help us understand the Fourier expansion at t = 0. Using this theorem, you will 
address question (5) in your homework. 

Theorem 4.1 (Some basic facts from Fourier analysis). If f(x) is a function such that 
I^Il 2 ([o i]) d ~ fo \f( x )\ 2 dx < °°> then f(x) can be Fourier- expanded as f(x) = £m=i A m sin(m7r.T), 
where A m = 2 f^ Q ^ f(x) sin(m7rx) dx. The infinite sum converges in the sense that 

N 

(4.0.14) I/- Y, ^mSin(m7rx)| L 2 ([0il]) -»■ as N -> oo. 

771=1 

We also have the Parseval identity 



(4.0.15) 



L2([0,l]) 



Y A^j sin(m7rx)|| 2 /r 0)1 -|s - Y T^^m 

777=1 A 



777=1 



Note that (4.0.15) is an "infinite dimensional Pythagorean theorem." 



Furthermore, if f is continuous on [0,1], then for any subinterval [a,b] c (0,1), 



(4.0.16) 



N 



\f - Y A rn sm(mirx)\ 



C°([a,6]) 



as N -*■ oo, 



777=1 



i.e., the convergence is uniform on any closed subinterval [a, 6] of the open interval (0,1). 



Exercise 4.0.1. Many extensions of Theorem |4.1| are possible. Read Appendix A of your textbook 
in order to learn about them. 
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